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Abstract 


The  structure  of  electrostatic 


oscillations  is  determined  for  a 


cold  plasma  with  a  aero  order  density  gradient  along  the  magnetic  field. 

A  set  of  well-behaved  structures,  bounded  at" plasma  resonance,  exists  for 
discrete  values  of  perpendicular  wavenumber.  Under  these  conditions  the 
plasma  acts  as  a  waveguide.  A  Green's  function  suitable  for  sources 


located  in  the  overdense  region  of  the  plasma  is  constructed  from  these 
waveguide-like  solutions . 
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I.  Introduction 


Magnetized  plasmas  with  a  zero  order  density  gradient  (7Nq)  along  the 
magnetic  field  (Bq)  can  support  a  rich  variety  of  electrostatic  oscillations. 

The  variety  of  mode  structures  in  this  geometry  arises  from  the  dose  juxta¬ 
position  of  the  plasma  (w^)  and  upper  hybrid  (co^)  resonances  for  wave  fre¬ 
quencies  (o  greater  than  the  electron  gyrofrequency  Q,  as  sketched  in  Fig.  1. 

The  majority  of  theoretical  studies  have  concentrated  on  two  extreme  cases. 

Either  the  plasma  is  unmagnetized  (Bq“0)  and  the  electrostatic  structure 
associated  with  the  plasma  resonance  is  studied,  or  the  density  gradient  is 
perpendicular  to  the  magnetic  field  (VN^B^O)  and  the  electrostatic  struc- 

tures  around  are  investigated.  Each  of  these  two  cases  has  received 
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considerable  theoretical  *  and  experimental  *  attention  and  numerous  interesting 
effects  have  been  found  for  each  case. 

In  the  geometry  considered  here  (VNQ  | |  Bq)  the  effects  of  both  reso¬ 
nances  must  be  considered  simultaneously.  In  fact,  WKB  analysis  of  a  wave 
propagating  obliquely  to  Bq  indicates  that  the  upper  hybrid  resonance  is 
actually  a  cutoff  for  parallel  wavenumber  (k„-0)  while  the  plasma  resonance 
is  a  true  resonance  (in  that  k„  ■+■  00 ) .  Due  to  the  coupling  between  the  m 

P 

and  <i>uh  points  and  the  need  to  consider  mode  conversion  near  the  point, 
the  general  problem  of  wave  excitation  by  an  external  source  is  rather  com¬ 
plicated  both  physically  and  mathematically.  Although  we  have  made  some 
progress  towards  a  solution  of  the  general  problem  we  concentrate  our  atten¬ 
tion  here  on  a  rather  specific  but  interesting  aspect  of  the  electrostatic  mode 
structures.  Specifically,  we  concentrate  on  a  description  of  new,  purely 
cold,  mode  structures  which  ere  supported  by  the  plasma  for  a  discrete  set 
of  perpendicular  wavenumbers  satisfying  the  condition. 


(1) 


kxAz  -[1  -  (0/fl))2]Js(2j  +  1) 
j  ■  0,1 ,2, « » • 

Where  As  refers  to  the  distance  between  the  upper  hybrid  and  plasma  resonance 
points  and  0  is  the  electron  gyrofrequeney.  These  nodes  are  trapped  in 
the  region  between  the  and  points.  The  wave  energy  propagates  across 

the  magnetic  field  and  the  region  behaves  as  a  waveguide.  We  demonstrate  that 
exact  solutions  to  the  cold  plasma  differential  equation  can  be  found  that 
are  well  behaved  when  the  condition  given  in  Eq.  (1)  is  satisfied.  With  these 
solutions  we  construct  a  Green's  function  suitable  for  sources  located  in 
the  overdenae  region  of  the  plasma  (u  <  u^) . 

Aside  from  the  intrinsic  interest  of  the  quantized  node  structures, 
their  existence  can  be  of  relevance  to  the  understanding  of  resonant  absorp¬ 
tion  of  radio  waves  in  the  auroral  ionosphere  or  in  certain  laboratory  devices. 
In  addition  this  study  gives  a  clear  example  of  a  node  structure  which  de¬ 
pends  explicitly  upon  the  plasma  lnhomogenlty  and  is  not  contained  in  the 
WKB  description  which  typically  forms  the  basis  for  physical  intuition. 


II.  Mode  Structure 


Consider  an  electrostatic  oscillation  propagating  In  an  Inhomogeneous 
plasma  in  which  the  zero  order  density  gradient  points  along  the  magnetic 
field  with  a  waveform  given  by 

♦  »  $(z)  exp[i(kjX  -  mt)]  +  c.c.  (2) 

The  direction  along  the  gradient  is  denoted  by  z  while  x  denotes  the  direction 
perpendicular  to  it*  as  illustrated  in  Fig.  1.  The  spatial  structure  of 
$ (z)  is  determined  by  Poisson's  equation 

*,<««»;♦>  -  kiVi*  -  S(z)  (3) 

where  the  perpendicular  and  parallel  components  of  the  dielectric  tensor  are 
ex  and  e„,  respectively.  In  a  cold  plasma  with  varying  electron  density 

ex  -  1  -  wp2(z)/(w2  -  «2)  ;  e„  -  1  -  wp2(z)/u>2  (4) 

2  2 

where  «p  (z)  -  4ire  NQ(z)/m  with  e,m  the  electron  charge  and  mass.  In  Eq.  (3), 
S(z)  represents  a  source  of  oscillations  at  frequency  u  which  may  consist, 
for  example,  of  externally  driven  grids  Immersed  in  the  plasma  or  of  charge 
oscillations  induced  on  the  density  gradient  by  a  remotely  launched  electro¬ 
magnetic  wave.  Regardless  of  the  particular  nature  of  the  source  the 
system  response  can  be  found  by  solving  Eq.  (3)  with  an  idealized  delta 
function  driver.  Mathematically,  this  is  equivalent  to  solving  for  the 
Green's  function.  Once  the  Green's  function  has  been  found  the  system  response 
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to  an  arbitrary  source  can  be  calculated. 

The  Green's  function  for  Eq.  (3)  can  be  constructed  from  linear  combin¬ 
ations  of  the  solutions  to  the  corresponding  homogeneous  equation  (i.e.  S(z)»0) 
which  satisfy  the  proper  boundary  conditions.  To  obtain  solutions  to  the 
homogeneous  equation  we  consider  m  density  profile  which  varies  linearly 
with  z  in  the  region  of  interest.  The  dielectric  components  then  take  the 
form 

e„  -  z/L  ;  ex  -  (z/L  -  Y2)/(l  -  Y2)  (5) 

where  Y  '■  fl/w  and  L  is  the  density  scale  length.  The  point  z  -  0  now 
corresponds  to  plasma  resonance  u  ■  The  overdense  region  of  the  plasma 

corresponds  to  z  <  0,  and  the  upper  hybrid  cutoff  is  located  at  z  ■  LY2. 

Since  we  assume  a  >  fl,  Y  <  1  and  we  introduce  the  dimensionless  quantities 


C  -  2|kjz/(l  -  Y2)1* 

(6a) 

and 

B  -  |kjL  Y2/2(l  -  Y2)** 

(6b) 

so  that  the  homogeneous  form  of  Eq.  (3)  becomes 

^•(C^)  -  (C/A  -  B)*  -  0 

(7) 

Using  the  transformation 


♦(C)  -  f(C)  «p(-C/2) 


(8) 


aawat _ >«**»*»»»* _ * 


Eq.  (7)  becomes 


C-^-r  f  +  (1  -  5)  4£  +  <e  -  *>  f  -  0 

as2  dK 


Eq.  (9)  has  the  form  of  Kummer's  differential  equation  which  has  two  independent 
solutions  usually  denoted  by  U(a,b,£)  and  M(a,b,£).  The  two  independent  para¬ 
meters  a  and  b  here  have  the  values  a  ■  *5  -  B  and  b  ■  1 .  Since  the  Wronskian 
of  U  and  M  is  zero  for  certain  values  of  a  and  b  we  use  a  different  com¬ 
bination  of  linearly  independent  solutions  to  construct  the  Green's  function. 

4 

In  the  notation  of  Slater  ,  the  solutions  we  use  are 


yfi  •  U(#,b,£) 


-  e5U(b-a,b,-£) 


The  Wronskian  of  y5  and  ^  is.  for  real  £ 


W(y5  ,j^)  -  £  empg  +  iff  (a  -  b)] 


which  is  not  zero  for  any  finite  value  of  a  or  b. 

The  Green's  function  is  a  solution  to  the  equation  (with  a  •  1) 


fjtC^G«.C0)]  -  (£/4  -  S)G(£,£o)  -  -4iroL6(£  -  £q)  (12) 


where  o  is  the  effective  charge  density  per  unit  area  of  the  source. 


_•  k 
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A  solution  to  Eq.  (12)  can  be  constructed  using  linear  combinations 
of  the  solutions  to  the  homogeneous  Eq.  (7)  on  either-  side  of  the 
source.  We  denote  these  combinations  by 


G  «  c  ’  exp(-£/2)y  +  c  exp(-?/2)y  ;  5  >  5 

>  11  5  12  7  c 


G  -  c  exp(-€/2)y  +  c  exp(-5/2)y  ;  £  <  E 

<21  5  22  7  o 


The  solutions  G>  and  G<  must  satisfy  the  following  conditions  to  form 
a  physically  acceptable  Green's  function: 

1)  They  must  he  continuous  at  £  *  £ 


«,«.>  -  «<«„> 


2)  The  derivative  at  ?  ■  5  must  satisfy 

o 


dG  dG  . 

/ _ >  _ <  \  -4ircrL 

(o5  -  ~ 


3)  Both  G>  and  G<  must  vanish  as  |g|->  »  and, 

4)  G>  and  G<  must  be  bounded  for  real  £. 

To  address  point  3  we  consider  the  asymptotic  behavior  of  the  function  U. 
For  |c|  ■+>  » 

U(a,b,E)  *  (0"a  (17) 


Using  Eq.  (17)  the  asymptotic  behavior  of  the  components  of  G>  and  G<  are 


exp<-e/2)y5  -  expC-S/aXC)0"5* 

and 

exp(-5/2)y  -►  exp(+€/2) (-5)”^ 

7 

where  we  have  used  the  values  of  a  and  b  appropriate  for  Eq.  (9). 
it  is  clear  that  in  order  for  G>  and  G<  to  vanish  as  |c|  ® 

c  a  c  m  0 
12  21 

Thus  the  Green's  function  must  be  constructed  from 
-  cne*p<-5/2)y5  ;  5  >  C0 

G<  "  C  22  e*P<-€/2)y 7  i  5  <  5o 


(18a) 


(18b) 


From  Eq.  (18) 


(19) 


(20a) 

(20b) 


In  considering  the  boundedness  oft  the  solutions  (point  4)  the  behavior 
of  the  solutions  at  plasma  resonance  (€  "  0)  is  crucial.  For  the  case  under 
consideration  (b  ■  1)  the  behavior  of  the  functions  U(a,l,£)  for  E  «  1  is 

U(a,l,0  =  -  In  C/r(a)  (21) 

• 

The  logarithmic  singularity  for  solutions  near  the  origin  when  r(a)  is  finite 
apparently  indicates  the  buildup  of  charge  near  the  resonant  layer.  In  this 

I 


case,  the  solutions  G>  and  G<  of  Eq.  (20)  can  not  be  used  to  construct  a 
physically  acceptable  Green's  function.  This  situation  Indicates  that  the 
cold  plasma  description  is  inadequate  to  describe  the  physical  processes 
occurring  near  the  resonance  layer.  Thermal  terms  must  then  be  retained  to 
adequately  describe  the  plasma.  However,  Eq.  (21)  suggests  the  possibility 
of  obtaining  physically  acceptable  solutions  when  r (a)  »  0,  i.e.  whenever 
a  is  zero  or  a  negative  integer.  Physically  this  implies  it  is  possible  to  obtain 
a  purely  cold  mode  which  vanishes  at  Infinity  and  is  bounded  at  plasma  reso¬ 
nance. 

The  parameter  a  can  indeed  have  zero  or  negative  integer  values  for 
certain  perpendicular  wavenumbers.  Thus  we  expect  bounded  solutions  to  exist 
whenever  a  >■  ^  -  £  «  -j  (j  ■  0,1,2,...).  Using  Eq.  (6b)  to  express  this  con¬ 
dition  in  terms  of  physical  variables,  kx  must  satisfy 

|kjL  -  (2j  +  1)(1  -  yVVy2  (22) 

For  these  discrete  or  ‘"umtlzed  values  of  k±  the  Kummer  function  U(a,l,£) 
becomes  a  polynomial,  namely, 

U(-j,l,C)  -  (-l)Jj!  L ,  <«  (23) 

where  Lj(£)  is  the  zero  order  Laguere  polynomial  which  is,  of  course, 

finite  at  the  origin.  Since  the  parameter  $  is  a  positive  quantity  the  solution 
G<  which  is  proportional  to  U(£+*s,l  ,-C)  is  not  bounded  at  the  origin. 

Thus  in  addition  to  satisfying  Eq.  (22)  we  must  also  require,  for  a  bounded 
Green's  function,  that  the  source  be  in  the  overdense  region  of  the  plasma. 
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(l.e.  5  <  0)  as  illustrated  in  Fig.  2.  Thus  the  solutions  and  vanish 

o  >  < 

at  infinity  and  are  everywhere  bounded  for  a  *  -j  and  CQ  <  0,  in  which  case 


G>  -  c  exp(-£/2) (-l)^j !  Lj(£) 


(24a) 


G>  -  c  exp(-C/2)  1^(0  ;  C  > 

and 

G  ■  c  exp(-C/2)y  (24b) 

<  22  7 

G  ■  c  exp(-?/2)U(j+l,l,-C)  ;  €  <  £ 

^  2  o 

The  constants  c  and  c  can  be  found  using  Eqs.  (15)  and  (16).  From 
1  2 

Eq.  (15)  we  find  * 

cz  -  cx  ecp(U0|)  Lj(-|Co|)/U(j+l,l,Uo|)  (25) 

The  derivatives  of  G>  and  G<  at  ?  ■  -CQ  *  |£q|  are 

G;  -  Cj  «p(|Co|/2)  [-ILj(-|C0|)  +Lj(-|50|)]  (26a) 

G^  -  c2  «xp(-|50|/2)[yU(l+j,l,|C0[)  +  U’ (j+1 ,1 , UJ )]  (26b) 


C2  "  “4liaL  eJ£pCU0|/2](j!)  LJ(-|5o() 


(27b 


Using  Eqs.  (27a, b)  the  Green's  function  appropriate  for  sources  on  the  over- 
dense  side  at  discrete  values  of  |kj  is 

G>  -  -(j!)  4,raL  U(j+1 ,1 , |cj )  LjU)  «cp[-(5  -  CQ)/2]  ;  (28a 

5  »  «„  *  -I«„l 

G<  -  -(J!)  4ircL  L^-UJ)  U(j+l,l,-5)  exp[(£  -  5^/2]  i  (28b! 
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III.  Discussion 

It  has  been  demonstrated  that  for 'Certain  discrete  values  of  perpen¬ 
dicular  wavenumber,  it  is  possible  to  construct  a  Green's  function  for  electro¬ 
static  modes  in  a  magnetized  cold  plasma  with  a  longitudinal  density  gradient. 
The  sources  must  be  located  in  the  over dense  region  of  the  plasma.  An  ex¬ 
ample  of  the  Green’s  function  for  j  >3  in  Eq.  (22)  is  shown  in  Fig.  3. 

The  Green's  function  is  oscillatory  in  the  region  between  plasma  resonance 
and  the  upper  hybrid  cutoff  point  beyond  which  it  decays  exponentially  as 
expected  from  WKB  analysis.  In  the  overdense  region  of  the  plasma  the  Green's 
function  decays  exponentially  away  from  the  source.  Viewing  the  Green's 
function  as  representing  physically  a  grid  limners ed  in  the  plasma  driven  by 
a  constant  current  source  of  amplitude  -i laa  it  is  clear  that  the  plasma 
response  is  purely  capacitive.  The  capacitive  impedance  of  the  plasma,  which 
can  be  found  from  Eqs.  (28a)  or  (28b),  yields  an  effective  plasma  capacitance 
per  unit  area,  C/A,  given  by 

(C/A)-1  -  (JI)  4irL  Lj(-|5o|)  U(j+1 ,1 , \ZQ\ )  (29) 

This  expression  illustrates  that  the  plasma  response  depends  upon  the  zero 
order  density  gradient  scale  length  and  the  details  of  the  electrostatic  mode 
structure. 

At  this  point  it  is  worth  noting  that  the  retention  of  terms  resulting 
from  changes  In  the  parallel  component  of  the  plasma  dielectric  is  crucial  to 
obtaining  these  results.  Neglecting  the  term  in  Eq.  (3)  proportional  to 
3  e„  leads  to  erroneous  results  in  at  least  two  ways.  First  the  behavior 

m 
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o£  the  electrostatic  modes  near  resonance  would  be  found  to  diverge  at 
nonquantlzed  values  of  |kA|  as  z”*  rather  than  the  correct  ln(z)  divergence. 
This  illustrates  that  near  resonance  the  plasma  can  not  be  treated  as  locally 
homogeneous.  Second,  the  differential  operator  would  be  non-adjoint  and 
any  subsequently  derived  Green's  function  would  be  incorrect  and  more  diffi¬ 
cult  to  obtain  than  the  correct  one. 

The  Green's  function  presented  here,  which  is  derived  from  the  cold 
plasma  equations,  is  valid  only  for  discrete  values  of  at  which  the 
quantity  F(a)  is  infinite.  At  other  values,  of  |kA|  the  solutions  to  the 
cold  plasma  differential  equation  have  a  logarithmic  singularity  at  plasma 
resonance.  Actually,  however,  the  amplitude  of  the  electrostatic  oscillations 
remains  finite  at  resonance  because  the  plasma  has  a  non-zero  temperature. 
Temperature  dependent  effects  lead  to  the  excitation  of  Bolm-Gross  modes  near 
plasma  resonance.  These  thermal  modes  carry  energy  down  the  density  gradient 
resonance  and  thus  Unit  the  amplitude.  Such  a  mode  conversion  pro¬ 
cess  is  described  mathematically  by  the  fourth  order  differential  equation  that 
arises  from  Eq.  (3)  if  thermal  terms  are  retained  in  e„.  The  investigation  of 
the  solutions  to  the  fourth  order  equation  is  the  subject  of  another  paper5  but 
for  completeness,  we  note  here  that  the  ratio  of  the  amplitude  of  the  hot 

generated  near  resonance  to  the  cold  mode  away  from  resonance  is  pro¬ 
portional  to  l/r(a).  Thus  the  amplitude  of  the  hot  mode  is  zero  at  the 
quantized  values  of  |kj,  consistent  with  the  results  presented  here. 

As  we  have  indicated  the  source  of  excitation  must  be  in  the  over- 
dense  region  of  the  plasma.  Even  at  the  quantized  values  of  |kj  purely 
cold  solutions  can  not  be  constructed  for  the  source  in  the  underdense  region 


because  of  the  inherent  singularity  at  resonance  in  the  solutions  11(04^,1 ,-£) . 
The  problem  is  quite  clearly  not  symmetric.  Obtaining  a  suitable  expression 
for  the  Green's  function  for  >  0  again  requires  the  retention  of  thermal 
effects.  The  solutions  to  the  fourth  order  differential  equation  remain 
bounded  at  the  resonance  because  of  production  of  hot  modes  through  mode 
conversion.  However,  finding  the  Green's  function  in  this  case  is  much  more 
complicated  because  four  modes  rather  than  two  can  be  excited  in  the  plasma. 

We  are  currently  investigating  this  rather  complex  problem. 

The  suppression  of  mode  conversion  at  discrete  values  of  [kx|  means 
that  the  plasma  acts  essentially  as  a  loss-free  waveguide.  The  wave  energy 
is  transported  perpendicular  to  the  field  lines  and  density  gradient.  The 
long  wavelength  cold  modes  are  not  subject  to  Landau  damping  whereas  the  short 
wavelength  Bobm-Gross  modes  created  through  mode  conversion  are.  The  experi¬ 
mental  investigation  of  the  waveguide-like  property  of  the  plasma  can  be 
carried  out  in  present  laboratory  devices  or  the  effect  may  be  observable  in 
absorption  experiments  conducted  in  the  auroral  ionosphere. 

From  the  limited  study  reported  here  it  is  evident  that  fundamental 
differences  exist  between  resonant  absorption  in  a  magnetized  and  unmagnetized 
plasma.  In  the  umagnetlzed  plasma  short  wavelength  Bohm-Gross  modes  are 
excited  by  direct  coupling  of  external  electromagnetic  energy  to  plasma  reso¬ 
nance.  In  the  magnetized  plasma  excitation  of  Bohm-Gross  modes  can  take 
place  indirectly  through  mode  conversion  from  cold  modes  and  is  suppressed  for 
certain  angles  of  propagation.  Thus  the  introduction  of  a  magnetic  field  in 
the  plasma  produces  not  just  a  change  in  the  angular  dependence  of  Bohm-Gross 
excitation  but  introduces  a  fundamentally  different  physical  mechanism  for 
the  production  of  these  modes. 
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Schematic  of  geometry  considered  in  this  study.  The  plasma 
density  gradient  VNq  points  along  the  magnetic  field  Bq.  With 
to  >  ft,  the  plasma  resonance  point  at  is  near  the  upper  hybrid 
point,  «uh. 

Schematic  indicating  the  relative  location  of  the  source,  plasma 
resonance,  and  upper  hybrid  cut-off,  To  obtain  bounded 

cold  solutions  the  driving  source  must  be  located  in  the  over- 
dense  region  (to^  >  to) .  Nq  Is  the  plasma  density. 

Example  of  spatial  dependence  of  a  cold  waveguide  Green's  function 
for  J  ■  3.  The  location  of  the  source,  plasma  resonance,  and 
upper  hybrid  cut-off,  are  indicated.  The  mode  is  finite 
everywhere  and  oscillatory  between  the  tu^  and  ta^  points.  <pg  is 
an  arbitrary  scaling  factor. 
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